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Present day results for few-nucleon bound state and scattering observ- 
ables based on modern high precision nuclear forces are briefly reviewed. 
While in relation to NN forces of that type three-nucleon (3N) forces are 
mandatory for binding energies and for quite a few 3N scattering observ- 
ables their effect is rather small in two-nucleon correlation functions as 
demonstrated for 3 He and 4 He. 

The old idea of the Coulomb sum rule as a way to extract the pp 
correlation function is reconsidered and the need for more accurate data 
is pointed out. It appears to be an ideal case to probe properties of the 
density operator and the ground state wave functions without disturbances 
of final state interactions (FSI). 

In the 3N system below the pion threshold FSI is well under control and 
therefore the exclusive process 3 He(e, e'NN) is also a very good test case 
for correlated nuclear wave functions and electromagnetic current opera- 
tors. One specific kinematics is emphasized, which can lead to insights into 
the correlated ground state wave functions with little disturbance of FSI. 
Finally exclusive photodisintegration of 3 He is regarded, which appears to 
be promising to identify 3N force effects. 

PACS numbers: 21.45+v, 21.10-k, 25.10+s, 25.20-x, 27.10+h 



1. Introduction 

The first and simplest description of nuclear physics is based on the 
non-relativistic Schrodinger equation. On this stage one requires that NN 
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interactions are well tuned to NN data up to the pion threshold and are 
accompanied by 3N forces, which at least guarantee the correct 3 H binding 
energy. For electromagnetic probes current operators consistent with the 
nuclear interactions are needed. To get insight into the dynamics of the 
nuclear systems, reliable solutions of such a Schrodinger equation are nec- 
essary. But extensions to this picture are possible. A more advanced and 
possibly necessary dynamical picture would include relativity for instance 
in the instant form of a Hamiltonian formalismfl]. Here we restrict our- 
selves to a strictly nonrelativistic treatment. We employ the present day 
perfectly well tuned NN forces CD-Bonn[2], AV18[3], Nijmegen I and II[4]. 
They are mostly of phenomenological and local nature, with the exception 
of CD-Bonn, which is a slightly modified one-boson-exchange potential and 
highly nonlocal. As 3N forces we choose the Tucson-Melbourne (TM) 2n- 
exchange [5] and the Urbana IX [6] models. In the first case the off-shell 
it — N amplitude is used in a low momentum expansion, in the second case 
that amplitude is based on an intermediate A-excitation and a phenomeno- 
logical short range part. Two strengths parameters in the Urbana IX 3N 
force are adjusted to the 3 H binding energy and nuclear matter density. In 
the TM 3N force we adjust the cut-off parameter A for a strong form- factor 
parametrization separately for each NN force partner to the 3 H binding 
energy [7]. The 3N system is solved rigorously in the Faddeev scheme for 
bound and scattering states; the 4N bound state is equally precisely evalu- 
ated using Yakubovsky equations [7]. A survey of typical and current results 
for binding energies and scattering observables is presented in Section II. 

Wave function properties in form of bound state two-body correlation 
functions are shown in Section III. One approach of connecting them to 
observables is the Coulomb sum rule. We briefly review that topic in Sec- 
tion IV and point to necessary improvements in experiment and theory in 
order to achieve clear and convincing results in the future. 

Another approach in investigating correlations are electron induced two- 
nucleon emissions, which we study in a specific kinematics for the target 
nucleus 3 He. It is shown in Section V that FSI appear to be unavoidable 
(at least below the pion threshold) but can possibly be reduced to an easily 
accessible and restricted one. This might enable a search for initial state 
correlations in a rather controlled manner. 

Photon induced two-nucleon emission on 3 He as well as the pd break-up 
appear to be very promising to see 3N force effects. This is illustrated in 
Section VI. Finally we end with a brief outlook. 
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Potentials 



3 H) £( 3 He) £'( 4 He) 



Nijm 93 
Nijm I 
Nijm II 
AV18 



7.668 
•7.741 
7.659 
7.628 
8.013 



7.014 
7.083 
7.008 
6.917 
7.288 



24.53 
24.98 
24.56 
24.28 
26.26 



CD-Bonn 



Exp. 



-8.482 



-7.718 



-28.30 



Table 1. 3 H, 3 He and 4 He binding energy predictions for several NN potential 
models compared to the experimental values. All energies are given in MeV. 



A first necessary test of the dynamical picture are three-and four-nucleon 
binding energies. For five modern, high precision NN forces our theoretical 
results are shown in Table 1. We have taken into account charge indepen- 
dence and charge symmetry breaking as well as the mass difference of the 
proton and the neutron. Further electromagnetic interactions are included 
and also the isospin T = | admixtures. We see the by now well known 
under binding against the experimental values. Table 2 collects the indi- 
vidually adjusted A-parameters of the TM 3N force and of a modified one 
(TM'), which violates chiral invariance less than TM, and the resulting 3N 
and 4N binding energies. We also show the AV18 plus Urbana IX results. 
We end up with the interesting result that the theoretical a-particle binding 
energies for the different NN and 3N force combinations are rather close to 
the experimental value. There is a slight over binding, which leaves little 
room for the action of 4N forces. This can be quantified in the following 
way. The average attraction due to NN forces is 24.92 MeV or 88 % from 
28.30 MeV. The average additional attraction due to 3N forces is 3.9 MeV 
or 14 % and finally the average over binding is 0.5 MeV, which is 2 % of 
28.30 MeV. If as a conjecture this would be attributed to a repulsive 4NF 
then this shows a nice hierarchy in the importance of two- to many-body 
forces. Of course this is a temporary statement and can be modified in the 
future if more will be known about strengths and properties of 3N forces. 
In any case such a hierarchy is in agreement with the expectations of chi- 
ral perturbation theory [9, 10]. An overview of 3N and 4N binding energies 
is shown in Fig. 1, which documents the strong correlation among them, 
known as Tjon line [7, 8]. 

The Argonne-Illinois-Los Alamos collaboration has explored with the 



2. Three-and Four-Nucleon Systems 
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Potentials 


A [m w ] 


E( 3 R) 


£( 3 He) 


£( 4 He) 


CD-Bonn+TM 


4.784 


-8.478 


-7.735 


-29.15 


AV18+TM 


5.156 


-8.478 


-7.733 


-28.84 


AV18+TM' 


4.756 


-8.448 


-7.706 


-28.36 


AV18+Urbana IX 




-8.484 


-7.739 


-28.50 


Exp. 




-8.482 


-7.718 


-28.30 



Table 2. 3 H, 3 He and 4 He binding energy predictions for several NN and 3N 
potential models compared to the experimental values. All energies are given in 
McV. 




Fig. 1. Tjon-line: a-particle binding energy predictions E( He) against the predic- 
tions for the 3 H binding energy for several interaction models. Predictions without 
(crosses) and with (diamonds) 3N forces are shown. The experimental point is 
marked by a star. The line represents a least square fit to NN force predictions 
only. 



help of the Greens-Function Monte Carlo method the low energy spectra of 
light nuclei up to A=8 [11]. We see in Fig. 2 the pure AV18 predictions, 
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Fig. 2. Spectra of light nuclei: experimental data (red, solid), AV18 (blue, dotted), 
AV18+Urbana IX (green,dashed). 



which are rather far away from the data and the impressive shift of theory 
towards the data by adding the Urbana IX 3N force. But we still observe 
deviations, which can be resolved by including additional terms in the 3NF 
[12] . Further important tests of the nuclear Hamiltonian without and with 
3N forces are scattering processes. For three nucleons solutions for the 
continuum are by far most developed. The Faddeev scheme [13] and the 
hyperspherical harmonic method [14] provide very accurate solutions. We 
illustrate the state of art with several cross sections and refer the reader 
to [13] and more recent papers [15, 16, 17] for a larger overview and for the 
very many spin observables, which probe our present day understanding 
of the dynamics in a very sensitive manner. Fig. 3 shows the nd total 
cross section, which below about 100 MeV is nearly perfectly described 
even without 3N forces. Only at the higher energies small discrepancies 
appear, which are however significantly reduced including 3N forces. At 
very low energies the inclusion of the pp Coulomb force is under control. For 
a survey on the beautiful agreement of the angular distribution in elastic 
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Fig. 3. The total nd cross section. Comparison of the data [33] with various inter- 
action models. 



pd scattering with precise data we refer to [14]. Thereby 3N force effects 
are tiny. This remains true also at somewhat higher energies as shown in 
Fig. 4, where the theory does not include the Coulomb forces nor 3N forces. 
At 65 and 135 MeV theory based on NN forces only clearly underestimates 
the data in the minima whereas the inclusion of 3N forces leads to a very 
good agreement. This is shown in Figs. 5 and 6. Thereby the results, as 
shown by two bands, are very stable under exchange of NN and 3N force 
combinations (which of course describe the 3 H binding energy). Finally we 
show in Figs. 7 and 8 some break-up cross sections along the kinematical 
locus as a function of a suitably defined arclength S [13]. In general the 
agreement is good and at those low energies 3N force effects studied up to 
now are insignificant. 
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A prerequisite to a theoretical analysis of 3 He(e, e'NN) data is the un- 
derstanding of the d(p,NN) reaction in the full phase space. Therefore 
47r-measurements of the latter process, not only at certain selected regions 
in phase space, are quite important to test the theory, before conclusions 
can be drawn from an analysis of the 3 He(e, e'NN) reaction. 

However, overall one can already say now that the dynamical picture 
with high precision NN forces and adjusted 3N forces works reasonably 
well (with room for improvements) and provides a good basis to analyze 
electromagnetically induced processes. 

3. Two-Nucleon Correlation Functions 

Based on fully converged solutions of the Faddeev-Yakubovsky equations 
we present two-nucleon correlation functions for ground state wave functions 
of the simplest type, namely averaged over two-body partial wave 

states: 
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E, =65 MeV 




Fig. 5. Differential cross section in elastic Nd scattering at 65 MeV. The light (dark) 
shaded bands are NN force only (NN+3NF) predictions for various interactions. 
pd data (crosses) are from [19] and nd data (open circles) are from [20]. 



ZJ + 1 M 

Here is the operator of a pair distance and Pij the projector on a pp or 
np pair. They are plotted (with arbitrary overall normalization) in Fig. 9 
for d, 3 He and 4 He and choosing two NN potentials AV18 and CD Bonn. 
At short distances below about 1 fm the predictions of the two potential 
are quite different. The curves are very much similar for the three nuclei, 
what suggests that this will essentially remain true also for heavier systems. 
They all peak at about 1 fm. Fig. 10 compares directly predictions for 
the three nuclei, all normalized to each other in the peak value. We see a 
nearly perfect overlap except in the tails, where the difference in separation 
energies show up. Finally the addition of 3N forces has no visible influence 
up to the radii somewhat larger than 1 fm, where binding effects have to 
appear. 
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Fig. 6. The same as in Fig. 5 for E= 135 MeV. pd data (crosses) are from [21] and 
(circles) from [22]. 



Since wave functions are no observables, in consequence the wave func- 
tion property C(r) is not an observable either. Wave functions enter into 
observable response functions as they occur for instance in electromagneti- 
cally induced processes, but they are accompanied by final state continuum 
wave functions, which are also correlated and, very importantly, they come 
together with current operators, which should be consistent to the nuclear 
forces. Consistent ingredients of the nuclear matrix elements should lead 
to the same observables (response functions) even if different NN force pa- 
rameterization have been chosen. At present mostly AV18 NN forces and 
consistent currents are being used. Thus work remains to be done to un- 
derstand possible model dependences. 
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10 20 

S [MeV] 




Fig. 7. Five-fold differ- 
ential Nd breakup cross 
section along the kine- 
matical locus including 
the quasi-free scattering 
condition, pd data at 19 
MeV are from ref. [23] 
and at 65 MeV from 
ref. [24]. 



4. The Coulomb Sum Rule 

It has been known for a long time [27, 28] that the Coulomb sum rule 
is an approach to isolate the Fourier transform of the correlation function 
C(r). Let's define the Coulomb sum as 

duR L (u,\Q\) (2) 

'mm 

Then using the standard expression for the longitudinal response func- 
tion Rl in inclusive electron scattering together with the closure relation 
one easily finds [29] 

Sl = \ ]T(*JM|ptp|*JM} - l -Y. \(*JM\p\*JM)\ 2 (3) 

M M 

This very nice intermediate result shows that the final state interactions 
have been totally removed and only ground state expectation values remain. 
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Fig. 8. The same as in Fig. 7 for two cases: two nucleons leave with equal momenta 
(up) and three nucleons leave with equal energies under pairwise angles of 120° 
(down), pd data at 13 MeV are from ref. [25], at 65 MeV from ref. [26]. 

Separating the density operator p into single particle and two- and more- 
particle parts one easily arrives at 

S L = ZG P E (Q) 2 + NG%(Qf - Z 2 Fl{Q) + C(Q) + C(Q) (4) 
where G p £ n are the p, n electric form factors ( neglecting the time component 
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of the four vector dependence in Q 2 ), is the elastic charge form factor 
of the target nucleus, C(Q) arises from the two- and more-particle densities 
and 
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Fig. 10. Comparison of C(r)'s for d, 3 He and 4 He (left) and for AV18 against 
AV18+TM (right). 



is the quantity related to C(r) suitably augmented by the nucleon electro- 
magnetic form factors [29]. 

There are certain obstacles to be overcome. The integral in Eq. 2 requires 
an extrapolation to catch all of the integrand above the quasi elastic peak. 
Up to now the available data leave too much room for ambiguities in the 
extrapolation [29] and precise data at some more higher w-values would 
be welcome. Also the first three terms on the right hand side of Eq. 4 
cancel strongly [29] leading to enhanced error bars for C and C. Therefore 
high accuracy measurements of Rl are needed. Above all this refers to the 
nuclei 3 He and 4 He, for which at present the most precise evaluation of the 
ground state wave functions are possible. On the theoretical side as has 
been shown in [30] the relativistic corrections in the density operator and 
two-body pieces therein play an important role and cannot be considered 
as a small perturbation. This is an interesting challenge for theory and 
experiment and deserves a renewed effort despite the intensive work in the 
past. 



Another approach towards correlations is via electron induced two-nucleon 
emission on nuclei, here on 3 He. The most ideal situation would be that one 
nucleon absorbs the photon and receives its full momentum and all three nu- 
cleons leave the nucleus without any final state interaction (FSI). Then the 




(5) 



5. Exclusive Electron Scattering on 



3 He 
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Fig. 11. The expression in Eq. 6 against the relative momentum p for a pp pair 
(left) and for a np pair (right). Both refer to NN forces only, whereas the figure in 
the middle includes in addition 3N forces. 



measurement of the momenta of the two spectator nucleons would display- 
directly their momentum distribution in 3 He. For this special kinematics 
the total spectator pair momentum has to be zero and those two nucleons 
leave back to back and show directly the relative momentum dependence 
within a pair of nucleons in the target nucleus 3 He. Let us number the 
nucleons such that the knocked out nucleon is number 1. Then one would 
probe directly the expression related to the 3 He wave function 
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PWIA 



PWIAS : 



RESCATT : 




Fig. 12. The matrix elements of Eq. 12 contributing to PWIA, PWIAS and to the 
infinite number of RESCATT processes. The first diagram in the group RESCATT 
is the expression " tGo " . 



£ £ |*(p,<f=0)| 2 (6) 

M mim2m3 

where p = \{k2 — fe) and q = §(&i — ^(A?2 + £3))- The fcj, i = 1, 2, 3 are the 
individual momenta. This quantity is plotted in Fig. 11 for the two spectator 
pairs pp and np and the two NN forces AV18 and CD-Bonn. We see model 
independence below about p = 1 fm _1 and significant model dependence 
at higher p-values. 3N force effects are, like for C(r), quite insignificant as 
shown in Fig. 11 for a np pair. 

Unfortunately reality is far from that ideal situation. FSI interferes very 
strongly. The eightfold differential cross section has the well known form 
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Fig. 13. The restrictions for relative NN momenta p as a function of photon mo- 
menta \Q\ choosing E^= 140 MeV. 



d 8 a 

crMott {vlRl + v t Rt + vttRtt + vtlRtl) P (7) 



dk'dk' <m 1( m2dS 

where the v's and p are kinematical quantities, and the response functions 
R are 

Rl = \N 1 2 

R T = |iVi| 2 + \N-i | 2 

R TT = 2Re(N 1 Nl 1 ) 

R TL = -2i?e(iVo(iV 1 + iV_!)*) (8) 

Here enter the nuclear matrix elements TVj, which are the spherical compo- 
nents of 

i^ = 3(*(-)|^(Q)|*s He > (9) 

with 

f(Q)=f(Q,l)+f(Q,23) + ... (10) 

Note that both 3N states in Eq. 9 are fully antisymmetrized. Using a 
Faddeev decomposition, can be written as 

iv^ = n£ wias + n£ escatt (11) 



L UND printed on February 8, 2008 



17 



where 

Kwias = 3(^(1 + ^(3)1*3^) 
Kescatt = 3(^-1(1 + P)Pn (12) 

and \U^) obeys a Faddeev-like integral equation [34] 

\U") = tG (l + P)f(Q)\*3He) + tPG \Un (13) 

The bra-vectors of Eq. 12 is given by the free momentum eigenstates 

\*%) = (1 - P 2 3)\p)\q) (14) 

antisymmetrized in the pair (23), t is the NN t-operator and Go the free 3N 
propagator. Finally P is the sum of a cyclic and anticyclic permutation of 3 
objects. It is instructive to display the physical content of Eq. 12 graphically 
in Fig. 12. There the different treatments of the final state, PWIA, PWIAS, 
"tGo" and "full" are explained. It is easy to show that in PWIA one finds 

R L = GI(Q)\^{P^=Q)\ 2 

Rt = 7^^(Q)i^|^(p,g = 0)| 2 (15) 

and Rtt = Rtl = 0. Thus Rl and Rt provide access to the same quantity 
J2 |^| 2 up to known factors. The most simple approximate FSI is the action 
of t within the spectator pair of nucleons, called " iGo" below in Figs. 14-17 
Since we are restricting ourselves to a nonrelativistic framework the total 3N 
cm. energy should be below the pion threshold at 140 MeV. Then for given 
\Q\ the p- values are restricted as shown in Fig. 13. For three |Q|-values, 400, 
500, 600 MeV/c, we compare in Fig. 14 the quantity R L /(F?) 2 for PWIA, 
PWIAS, " tGo", and using the complete FSI (called Full), as a function of 
the relative momentum of an outgoing np pair. In PWIA pis related to the 
spectator pair (the above mentioned pair 23). Consequently it is assumed 
that the knocked-out nucleon is a proton. Of course for PWIA there is no 
|Q|-dependence and we see directly the expression given in Eq. 6. FSI has a 
tremendous effect for all three |Q|-values, however the simple approximation 
" tGo " improves going to the higher |Q|-values. Even for higher \Q\ values 
FSI has a tremendous effect on Rl (a reduction factor 10-100). At least for 
the higher \Q\ values the tGo approximation might be reasonable, but still 
one can see effects of the neglected FSI of spectator nucleons and hit proton. 
In case of Rt we show in Fig. 14 the ratio {2wi 2 n Rt) / {G P M Q) 2 . In addition 
some curves include it and p exchange contributions consistent to AV18. 
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Their effect is totally negligible. Also in this case the approximation " iGo" 
is a quite good representation of the full FSI. Since that FSI correction of 
first order in t is very well under control given the fact that the on-shell t 
is fitted to the NN data, an analysis of future data appears to be a rather 
reliable approach towards the quantity given in Eq. 6. There is nice stability 
under exchange of the NN forces AV18 against CD Bonn, as shown in Fig. 15 
for the example \Q\= 400 MeV/c. 

The situation is quite different if the hit nucleon is a neutron. We 
show in Fig. 16 the quantity Rl/{G%) 2 . In contrast to the case before, 
the approximate FSI treatment " tGo " is now totally different from the 
complete FSI result and thus would be totally misleading. This is of course 
due to the smallness of G%. For Rt, however, the " tGo " approximation 
and " Full " are not far from each other and moreover the |Q|-dependence 
is relatively weak. Nevertheless for both R 7 s the predictions are nearly 
independent from the specific choice of the NN force as shown in Fig. 17 
and therefore the model dependence is weak. 

All the curves in Figs. 14 and 16 refer to AV18 and the angle between p 
and Q has arbitrarily been fixed at 90°. At the other angles the relation be- 
tween the "full" result and the " tGo " approximation changes, but remains 
within the same order of magnitude. 

Based on these results one has to state that there is no way to access di- 
rectly in that low energy regime \ ^\ 2 - Nevertheless precise measurements 
would be extremely informative to test the whole dynamical picture, forces 
and currents. Choosing other kinematical conditions even in PWIA the 3 He 
bound state *S>(p, q) is probed in such a manner that both p and q vary. Un- 
der the prerequisite that the pd break-up process has been tested carefully 
against theory the 3 He(e, e'NN) reaction for general kinematics is a perfect 
tool to probe the remaining unknown ingredient, the current operator [31]. 



6. 3N Force Effects in Photo-Induced Disintegration of 3 He 

3N forces are required for a correct description of binding energies. Will 
they also play a role in electromagnetically induced processes ? We already 
started a first investigation for pd capture processes in [32] and would like 
to show new results for photodisintegration of 3 He, specifically at higher 
energies than the one considered in [32]. To that aim the Faddeev-like 
integral equation in [34] has to be modified. The nuclear matrix element for 
3 He (7, NN) 

N» = (^Z>\f\* 3He ) (16) 
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can be written as 

N» = ±{<S>%\(l+tG )P\U») (17) 
where \U^) obeys the Faddeev-like integral equation 

\U») = (1 + P)f\**He) + (tG P + + P)V} 1] G (1 + tG )P)\un (18) 

We use the Siegert theorem as described in [32]. This includes some of the 

exchange currents. The operator vj 1 ^ is that part of the 3N force, which is 
symmetric under the exchange of particle 2 and 3. We scanned the whole 
phase space comparing the full break-up cross section evaluated with and 
without 3N force. As an example we display in Fig. 18 the 61,62 regions 
where 3N force effects are larger than 40 %. The relative azimuthal angle 
is < 50°. 

As an illustration we show out of that phase space region three more 
or less arbitrarily selected break-up cross sections along the S-curve in 
Figs. 19-20. The peak in the middle is caused by small relative momenta in 
one pair (FSI peak). In the other peaks the 3 He wave function is probed at 
small momenta. Measurements should validate or invalidate these sort of 
predictions. Also in the pd break-up process 3NF effects are clearly visible 
as shown in Fig. 21 for various photon energies. 

7. Outlook 

One main theoretical challenge is to establish an electromagnetic cur- 
rent operator, which is consistent to nuclear forces. Only then response 
functions for electromagnetically induced processes can be put on a firm 
ground. Precise experimental data on two-nucleon emissions induced by 
real and virtual photons on 3 He will be a very important test ground to 
probe nuclear forces, correlated wave functions and currents. At present 
the 3N system is the only case where FSI is fully under control (below the 
pion threshold) and appears therefore especially promising to probe our 
present day understanding of nuclear dynamics. 
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Fig. 14. R L /(F?) 2 (up) and (2m 2 N R T )/ (G P M Q) 2 (down) against the np relative 
momentum p for various treatments of the final 3N state and the |Q|-values 400, 
500 and 600 MeV/c. MEC effects are negligible as shown. 
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Fig. 15. R L /(F?) 2 (up) and (2m 2 N R T ) / (G P M Q) 2 (down) against the np relative 
momentum p for various treatments of the final 3N state and two different potentials 
at one \Q\- value. 
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Fig. 16. R L /(G n E ) 2 (up) and (2m 2 N R T ) / \G n M Q) 2 (down) against the pp relative 
momentum p for various treatments of the final 3N state and some |Q|-values. 
MEC effects are negligible as shown. 
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Fig. 17. R L /(G n E ) 2 (up) and (2m%R T ) / \G n M Q) 2 (down) against the pp relative 
momentum p for various treatments of the final 3N state and two different potentials 
at one \Q\- value. 
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Fig. 18. The 6>i-02 re- 
gions in 3N phase space, 
where 3N force effects 
are larger than 40 %. 
The relative azimuthal 
angle is </>i 2 < 50°. 
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Fig. 19. Three-body differential photo-disintegration cross sections 

a? a / (dVLid£l 2 dS) [lm 2 /(sr 2 MeV)] of 3 He along the kinematical locus for dif- 
ferent combinations of angles. 
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Fig. 20. The same as in Fig. 19 for a specific breakup configuration, where two 
nucleons leave with equal momenta leading to a FSI peak. 
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Fig. 21. The pd photodisintegration breakup cross section da/dttd [ fm 2 /sr ] of 
3 He based on AV18 alone (dashed curve) and AV18+Urbana IX (solid curve) for 
various photon energies q in MeV. 
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